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ABSTRACT 

DM Cyg, a fundamental mode RRab star was observed in the 2007 and 2008 
seasons in the frame of the Konkoly Blazhko Survey. Very small amplitude light curve 
modulation was detected with 10.57 d modulation period. The maximum brightness 
and phase variations do not exceed 0.07 mag and 7 min, respectively. In spite of the 
very small amplitude of the modulation, beside the frequency triplets characterizing 
the Fourier spectrum of the light curve two quintuplet components were also identified. 
The accuracy and the good phase coverage of our observations made it possible to 
analyse the light curves at different phases of the modulation separately. Utilizing the 
IP method (Sodor. Jurcsik and Szeidl, 2009) we could detect very small systematic 
changes in the global mean physical parameters of DM Cyg during its Blazhko cycle. 
The detected changes are similar to what we have already found for a large modulation 
amplitude Blazhko variable MW Lyrae. The amplitudes of the detected changes in the 
physical parameters of DM Cyg are only about 10% of that what have been found in 
MW Lyr. This is in accordance with its small modulation amplitude being about one 
tenth of the modulation amplitude of MW Lyr. 

The pulsation period of DM Cyg has been increasing by a rate of fj = 0.091dMyr _1 
during the hundred-year time base of the observations. Konkoly archive photographic 
observations indicate that when the pulsation period of the variable was shorter by 
Appuls = 5 • 10~ 6 d the modulation period was longer by Ap nlo d = 0.066 d than today. 

Key words: stars: variables: other - stars: horizontal branch - stars: individual: 
DM Cyg - techniques: photometric - methods: data analysis - 



1 INTRODUCTION 

Utilizing our full access to an automatic 60 cm telescope we 
have obtained extended multicolour observations of many 
fundamental mode RR Lyrae variables showing light curve 
modulation (the Blazhko effect) during the past five years. 
Detailed analyses of some of our targets were already pub- 
lished in I Jurcsik et~all (|2005l , 120061 , 120081 . 120091 ). Our obser- 
vations are the first multicolour photometric data which are 
accurate and dense enough to allow not only to determine 



* Based on observations collected with the automatic 60 cm tele- 
scope of Konkoly Observatory, Budapest, Svabhegy 



the frequencies appearing in the Fourier spectra of the light 
curves, but also to find changes in the global mean physical 
parameters of the stars (L,T e ff,R) in different phases of 
the modulation. In order to extract changes in the phys- 
ical parameters exclusively from multicolour photometric 
data we devel oped an inverse photometric B aade-Wesselink 
method (IPM: ISodor. Jurcsik fc Szeidlll2009l ). The IPM was 
successfully applied to determine changes in the mean global 
physical parameters of MW Lyr during its Blazhko cycle in 
lJurcsik et all (|2009h . 

DM Cyg, a short period (P = 0.42 d) RRab star was 
selected to be observed in the frame of the Konkoly Blazhko 
Survey because it had been announced to show phase modu- 
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Table 1. Konkoly photometric observations of DM Cyg. 



HJD 


mag* 


detector 


filter 


2400000+ 








54265.45479 


-0.249 


CCD 


V 


54265.46045 


-0.404 


CCD 


V 


54265.46611 


-0.501 


CCD 


V 


* Relative magnitudes for the CCD and photoelectric data. 



Table 2. Normal maximum timings derived from the Konkoly 
photographic, photoelectric and CCD observations. 



.ID 



Normal maximum timings obs. 



2427667 - 2429226 
2434186 - 2436373 
2443747 - 2444133 
2454265 - 2454407 
2454597 - 2454711 



2429078.482 

2435342.347 

2443777.333 

2454334.3467 

2454661.4210 



Pg 

Pg 

pe 
CCD 
CCD 



lation with a period of 26 d by Lvsova fc Firmanvukl (|l980t ). 
but neither the NSVS data jWozniak et all bOO^ Tnor the 
complete data set of its maximum timings show phase mod- 
ulation with this period (|S6dor fc Jurcsikll2005l ). 

Though DM Cyg is a relatively bright (V = 10-11 mag) 
RR Lyrae star, complete, accurate, multicolour light curve 
of its pulsation has never been published. Its period change 
was, however, regularly monitored by different groups of ob- 
servers. The GEOS databastf] lists 260 maximum times of 
DM Cyg between 1900 and 2008. 

In the present paper we publish our extended CCD ob- 
servations of DM Cyg and the results of the analysis of its 
light curve modulation. Archive photographic and photo- 
electric Konkoly data are also processed. 



2 DATA 

CCD observations were obtained with the automated 60 
cm telescope of the Konkoly Observatory, Svabhegy, Bu- 
dapest equipped with a Wright Instruments 750 x 1100 
CCD camera and BVIc filters. Measurements were taken 
on 81 nights between July 2007 and Sept 2008. About 
3100 data points in each band were gathered. Exposition 
times were 200, 60 and 40 sec or a bit longer depend- 
ing on the sky transparency in the BVIc bands, respec- 
tively. Data reduction was performed using standard IRAF0 
packages. Aperture photometry of DM Cyg (21:21:11.548 
+32:11:28.71) and several neighbouring stars were carried 
out in order to check the stability of the photometry and 
the constancy of the comparison stars. The relative mag- 
nitudes of DM Cyg measured to the mean magnitudes of 
C1=GSC2.2 N0330220980 (21:21:30.731 +32:13.05.78), and 



1 http : //dbrr . ast . obs-mip . f r /maxRR . html 

2 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Figure 1. Delta V magnitudes versus Julian Date, data phased 
with the 0.419863 d pulsation and the 10.57 d modulation periods 
are shown in panels A, B and C, respectively. 



C2=GSC2.2 N03302207371 (21:20:59.483 +32:13.01.63) are 
used in the analysis. The Tycho B and V magnitudes of the 
comparison stars CI and C2 are B T = 12.705, V T = 11.799, 
and B T = 13.081, V T = 12.150, respectively (|Hog et all 
|2000| ). Second order extinction correction of the data were 
applied in the B band. Relative magnitudes are transformed 
to standard BVIc magnit udes. More details o n the reduc- 
tion procedure are given in I Jurcsik et al.l (|2008l ). 

Archive photoelectric and photographic data obtained 
with the 60 cm telescope in 1978 and with a 16 inch astro- 
graph between 1934 and 1958 are also utilized. Photoelectric 
observations were obtained on 4 nights, the BV magnitudes 
were measured relative to GSC2.2 N03302207371. The pho- 
tographic measurements comprise data from 40 nights, the 
plates were evaluated using GSC2.2 B magnitudes of the 
surrounding stars. 

Table [1] shows a sample of the Konkoly photometric 
data of DM Cyg. All the photometric observations are avail- 
able in the online version of the journal as Supplementary 
Material. 

Seasonal normal maximum timings from each data set 
have been determined for five epochs, these data are listed 
in Table [2] 



3 RESULTS 

The accurate CCD observations reveal that the light curve of 
DM Cyg is not stable. It shows small amplitude modulation, 
the variations in maximum brightness and phase are only 
about 0.07 mag and 0.005 d (7 min, 0.012 pulsation phase), 
respectively. 

The V light curve of DM Cyg and the results obtained 
folding the data with the pulsation and modulation periods 
are shown in the three panels of Fig. [T] The phased resid- 
ual light curve after the removal of the mean pulsation light 
variation is given in Fig. [2] This figure shows the differences 
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Figure 3. Maximum V brightness and maximum phase values phased according to the Blazhko period are plotted in the top-left and 
bottom-left panels, respectively. Their typical errors are 0.005 mag and 0.002-0.004 phase. The right-hand panel shows the maximum 
brightness vs maximum phase data. The progression of data during the Blazhko cycle in this diagram is anti-clockwise. Dashed lines 
connect the mean values in 10 phase bins of the modulation. 
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Figure 2. Residual V light curve of DM Cyg after removing the 
pulsation components from the data according to the light curve 
solution given in Table \3\ Residuals of the smallest (stars) and 
the largest (diamonds) amplitude phases are shown with large 
symbols in the left-hand and right-hand panels, respectively. Note 
that the residuals have larger amplitude at around the middle of 
the rising branch of the pulsation (phase = —0.1) than around 
pulsation maxima (phase = 0) due to phase modulation. 

of the amplitude of the modulation at different phases of the 
pulsation. The well observable amplitude of the modulation 
is concentrated to a narrow phase range of the pulsation 
around minimum, rising branch and maximum. A similar 
behaviour has been already found in RR Gem and SS Cnc, 
in two other small modula tion amplitude Blazhko variables 
jjurcsik et al.ll2005l , [200(3 ). The amplitude of the the resid- 



ual variation is the largest at around the middle of the rising 
branch of the pulsation light curve (phase = —0.1) indicat- 
ing phase modulation of the rising branch as well. There is a 
difference between the phase of the amplitude modulation of 
the maximum of the light curve and the phase of the phase 
modulation of the same feature. In particular the maximum 
positive displacement of the timings of the maximum bright- 
ness precedes the occurrence of the brightest maximum by 
about 0.25 Blazhko phase as it is shown Fig. [3] 

The elements of the pulsation and the modulation are: 

T maxpuls = 2 454 312.514 [HJD] + 0.419863 ■ E puls , (1) 

and 

TmaxBi = 2 454 312.514 [HJD] + 10.57 ■ E m . (2) 

-Epuls and Ebi denote the epoch number of the pulsation and 
modulation cycles, respectively. 

The pulsation and modulation periods correspond to 
the mean values of the frequencies of the light curve so- 
lutions of the B, V and Ic data using 'locked' frequency 
solutions of the kfo pulsation and fc/o ± / m , and / m modu- 
lation frequency components (i.e., the modulation side lobe 
components are at the positions of the linear combination 
frequencies) . 

Data analysis was performed u sing the differ ent appli- 
cations of the MUFRAN package l |Kollath|[i99ch . a linear 
combination fitting program developed by A. Sodor, and 
the linear and nonlinear curve fitting abilities of gnuplo10. 

3 http://www.gnuplot.info/ 
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Figure 4. The left-hand panels show the spectral window, the amplitude spectrum and the residual spectrum after the removal of the 
pulsation components for the CCD V data. The middle and right-hand panels show residual spectra in the vicinity of the 10/o,H/o, 
and 12/o (fc = 10, 11, 12) and of the 25/o, 26/o, and 27/o (k = 25, 26, 27) frequencies, respectively. The spectrum has been prewhitened 
for the pulsation components up to the 24th order, and in the right hand panels for the detected modulation components, too. These 
spectra indicate that modulation frequency components are present only in the first 10-11 orders, while the pulsation can be accurately 
described if harmonic frequency components are taken into account up to as high as the 27th order. For clarity, vertical grids denote the 
positions of the fc/o, fc/o + fm and kfo — f m frequencies. 
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Figure 5. Fourier spectra of the V and B residual data in 
the vicinity of /o and 2/rj. The pulsation and the fc/o ± f m 
modulation frequencies have been removed. The /o — 2/ m and 
2/o — 2fm quintuplet components appear in these residual spec- 
tra with ~ 0.001 mag amplitude. 
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3.1 The light curve solution 

The Fourier amplitudes and phases of the pulsation and 
modulation frequency components identified in the spectra 
of the BVIc light curves of DM Cyg are summarized in 
Table [3] The Fourier decompositions use sin terms and the 
initial epoch corresponds to one of the brightest maxima of 
the modulated light curve, TO = 2454312.514. The errors of 
the amplitudes are ~ 0.0002 mag, the errors of the phases for 
the time transformed V data (see the details later) are given 
in the last column. The pulsation components are detectable 
up to the 27th order, while the modulation side frequencies 
(fc/o±/m) diminish at around the 11th order as documented 
in Fig. [4] The modulation frequency f m is present without 
question in each of the B,V,I spectra with similar ampli- 



Figure 6. Amplitude decrease of the pulsation and modulation 
side lobe frequencies at different orders shown as amplitudes nor- 
malised to the amplitudes of the first order (k=l) components. 
The decrease of the amplitudes of the harmonic components of 
the pulsation is exponential like, while the decrease of the ampli- 
tudes of the side lobe frequencies is more linear, with irregular 
character in the low orders. 



tudes as the 4th order negative and the 9th order positive 
modulation side lobe components have. Quintuplet frequen- 
cies are detected at fo — 2f m and 2/o — 2/ m as shown in 
Fig. El 

The 51 frequencies listed in Table[3]fit the data with 6-8 
mmag residual scatter in each band, which is about the level 



Photometric study of the Blazhko star DM Cyg 5 



Table 3. Fourier amplitudes and phases of the pulsation and modulation frequencies in DM Cyg. Initial epoch corresponds to maximum 
of both the pulsation and the modulation light variations as given in Eqs. 1 and 2. 



frequency B V I c V t * t a v{V tt ) 









amp 




60 


amp 




60 
r 


amp 




60 
T" 


amp 




60 
r 






fm 




0.094607 


0.0026 


1 


.678 





0017 


1 


.637 





,0014 


1 


,517 


0.0017 


1 


,644 





,107 


fo- 




2.192515 


0.0015 


3 


363 





.0008 


2 


.989 





,0005 


3 


,939 


0.0008 


2 


,697 





,211 


fo- 


./ m 


2.287122 


0.0036 


2 


.694 





0022 


2 


618 





0011 


2 


839 


0.0048 


3 


,675 





,037 


fo 




2.381729 


0.4903 


3 


.944 





.3532 


3 


.874 





2170 


3 


,670 


0.3532 


3 


,874 





001 


fo + 




2.476337 


0.0132 


3 


.822 





.0096 


3 


.854 





,0061 


3 


,975 


0.0057 


1 


,043 





,032 


2/o - 


" 2/ttt, 


4.574244 


0.0014 





.664 





.0008 





.275 





,0004 





,976 


0.0008 





,573 





,219 


2/o - 


" /m 


4.668851 


0.0035 


2 


.355 





0026 


2 


.404 





,0012 


2 


513 


0.0043 


3 


757 





,042 


2/o 




4.763459 


0.2598 


1 


.109 





1919 


1 


.090 





,1189 


1 


,029 


0.1920 


1 


,090 





,001 


2/oH 


~ /m 


4.858066 


0.0105 


3 


.511 





.0080 


3 


.572 





,0049 


3 


,619 


0.0038 


3 


,284 





,047 


3/o ^ 


/m 


7.050581 


0.0050 


3 


173 





.0034 


3 


187 





0024 


3 


,270 


0.0053 


1 


,058 





,034 


3/o 




7.145188 


0.1545 


1 


475 





1160 


1 


.476 





,0728 


1 


,466 


0.1162 


1 


,476 





,002 


3/o H 


- fm 


7.239795 


0.0100 


3 


928 





.0079 


3 


.988 





0048 


4, 


,003 


0.0040 


3 


772 





,044 


4/o- 


~ fm 


9.432310 


0.0024 


3 


.179 





.0019 


3 


.237 





,0012 


3 


,230 


0.0027 


1 


,562 





,065 


4/o 




9.526917 


0.0816 


1 


.970 





.0628 


4. 


.971 





,0396 


1 


,972 


0.0629 


1 


,970 





,003 


4/o H 


~ fm 


9.621525 


0.0086 


1 


.464 





.0063 


1 


.437 





0042 


4, 


466 


0.0036 


1 


222 





,050 


5/o- 


~ fm 


11.814039 


0.0019 


3 


.800 





.0016 


3 


.799 





,0010 


3 


,646 


0.0027 


1 


,826 





,067 


5/o 




11.908647 


0.0510 


5 


.142 





.0391 


5 


.171 





,0247 


5 


,204 


0.0392 


5 


,173 





,004 


5/oH 


~ fm 


12.003254 


0.0069 


1 


868 





.0053 


1 


889 





,0035 


1 


,862 


0.0030 


4, 


,874 





,060 


6/o - 


~ fm 


14.195768 


0.0012 


4, 


312 





0011 


1 


130 





,0008 


1 


,502 


0.0019 


5 


,325 





,094 
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5 


.592 





.0255 


5 


.617 
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5 
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5 
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,007 


6/oH 
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0.0056 


5 


.159 





.0045 


5 


.169 





,0027 


5 


102 


0.0028 


5 


,047 





,064 


7/o- 


~ fm 


16.577498 


0.0015 


5 


.147 





.0010 


1 


.585 





,0010 


1 


,748 


0.0015 


5 


,503 





,117 


7/o 




16.672105 


0.0185 


5 


.861 





.0146 


5 


921 





0088 


5 


,979 


0.0147 


5 


,925 





012 


7/o H 


~ fm 


16.766713 


0.0040 


5 


.601 





.0034 


5 


591 





,0020 


5 


,667 


0.0022 


5 


527 





,080 


8/o- 


~ fm 


18.959227 


0.0003 





.174 





.0006 


5 


.416 





,0004 


6 


,238 


0.0012 


6 


113 





,143 


8/o 




19.053834 


0.0114 


6 


253 





.0089 





.030 





,0050 





,167 


0.0089 





033 





,018 


8/o H 


~ fm 


19.148442 


0.0030 


5 


913 





0026 


5 


.954 





,0016 


6 


,054 


0.0018 


5 


,879 





,099 


9/o - 


~ fm 


21.340956 


0.0009 





,444 





0006 





.525 





0009 





135 


0.0010 





,647 





,172 


9/o 




21.435564 


0.0048 





.309 





.0037 





.482 





,0024 





,769 


0.0038 





,480 





,046 


9/o H 


~ fm 


21.530171 


0.0021 





.204 





.0019 





034 





0013 





,248 


0.0015 


6 


,239 





,113 


10/o 


~ fm 


23.722686 


0.0004 


1 


.793 





.0005 


1 


.018 





,0008 





,687 


0.0006 


1 


,077 





,289 


10/o 




23.817293 


0.0016 





.291 





.0012 





.760 





,0007 


2 


,221 


0.0012 





,780 





,141 


10/o 


H~ f m- 


23.911900 


0.0015 





.576 





0013 





.624 





,0009 





,652 


0.0011 





,573 





162 


H/o 




26.199022 


0.0006 


4. 


.556 





.0009 


1 


.183 





,0009 


3 


,499 


0.0009 


1 


,147 





192 


H/o 


+ fm 


26.293630 


0.0009 


1 


.438 





.0010 





.835 





0008 


1 


,374 


0.0011 





,818 





162 


12/o 




28.580752 


0.0019 


4, 


.559 





0014 


1 


.526 





0015 


1 


,265 


0.0015 


1 


462 





,116 
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4 
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1 
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1 
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1 
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5 
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5 
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5 
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5 
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5 
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6 


.073 
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5 


.883 
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5 
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5 


,845 
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6 
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6 
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6 
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6 
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,088 
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.187 





.0017 
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,062 


0.0017 





106 





,096 


19/o 




45.252857 


0.0016 





.452 





.0014 





,546 





,0015 





,444 


0.0015 
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111 


20/o 




47.634586 


0.0018 





.765 





.0014 





920 





0011 





,755 


0.0014 





,881 





,121 
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50.016315 


0.0009 
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.0012 


1 


,023 





,0011 


1 


,094 
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,947 





,132 
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.161 
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1 
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,0010 
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,384 
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1 
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23/ 
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.0010 


1 
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,658 


0.0010 


1 
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,177 
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57.161503 


0.0009 


2 


.064 
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2 


,306 





0006 


2 


,073 
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2 


,259 





,177 
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59.543233 


0.0006 


2 
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of the observational noise. No further periodic or stochastic 
variation in the residuals is detected. 

The decrease of the amplitudes of the pulsation compo- 
nents with increasing order is smooth and exponential-like. 
The amplitude decrease of the modulation components is 
however, different. The amplitudes of the low order mod- 



ulation components behave irregularly, while they show a 
linear decrease at higher orders (see Fig. IfoJ). Similar be- 
haviour of the amplitude decrease of the detected frequen- 
cies were demonstrated for RR Gem, SS Cnc a nd SS For 
(IJurcsik et al. I l2005l . 12006] ; iKolenberg et al.ll200Sl ) . 

Fig.[7]shows the V light curves belonging to 10 different 
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Figure 7. AV light curves in 10 phase bins of the modulation. All the observations are also shown with gray colour in each plot. Empty 
circles denote artificial data which are used to stabilize the Fourier fits if there is some gap in the observations. 



phase bins of the modulation. The first five Fourier param- 
eters of the 15th order harmonic fits to these individual V 
light curves are plotted in Fig. [8] Comparing the variations 
of the Fourier parameters of the l ight curves of DM Cyg 
with that of MW Lyr (Fig. 12 in I Jurcsik et al.l l2008h the 
following differences are conspicuous. The amplitudes of the 
detected changes in the Fourier parameters are about one 
order of magnitude smaller than the amplitudes of the light 
curve variations in MW Lyr. The amplitudes of the varia- 
tion of the amplitudes of the fo, ...5/o components decrease 
more drastically with increasing order in MW Lyr than in 
DM Cyg. This is most probably connected to that the pul- 
sation components are detected only up to the 12th order in 
MW Lyr, while they can be observed up to the 27th order in 
DM Cyg. Nevertheless the phase relations between the am- 
plitude modulation and the variation in the phase of the fo 
pulsation component are the same for the two stars, the de- 
tected changes in the epoch independent phase differences 
[v(/fei)] are, however, significantly different. While in DM 
Cyg the tp(fki) phase differences show sinusoidal variations 
with 0.25 (it/2) phase shift and with larger amplitude than 
<p(fo), in MW Lyr the amplitudes of the variations of the 
P(fki) components are smaller than the amplitude of <p(fo) 
and they do not show strictly regular behaviour with the 
modulation period. 

In lJurcsik et al l i|200Sl ) we have introduced a new 
method of analysing Blazhko variables' light curves. It was 
shown that if the times of the observations are corrected 
according to an appropriate time transformation defined by 
the variation of the phase of fo during the Blazhko cycle, 
then the modulation can be separated into phase and am- 
plitude modulation components. The time transformed data 
showed pure amplitude modulation for MW Lyr. Though the 
modulation amplitude of DM Cyg is very small, we have 
checked how a similar time transformation influences its 
modulation properties. Transforming the times of the ob- 



servations according to a continuous harmonic function of 
the tp(fo) data, the V observations folded with the pulsation 
period is compared to the folded light curve of the original 
data of DM Cyg in Fig. Enlarged plots of the middle of 
the rising branch are inserted in the figures. Though the dif- 
ference between the original and the time transformed data 
is small, the reduction of the phase modulation component 
in the time transformed data is indicated by the narrowness 
of the rising branch in the right-hand panel of Fig. [5] 

In the last three columns of Table [3] the Fourier ampli- 
tudes and phases and the errors of the phases are listed for 
the time transformed V data. The amplitudes of the mod- 
ulation side frequency components kfo + f m a,nd kfo — f m 
are about the same and the phases of the triplets in the dif- 
ferent orders show definite phase coherency according to the 
Fourier solution of the time transformed data (see Fig. 1 10p . 
As the phases of the components of the triplets are epoch de- 
pendent, this phase coherency holds only if the initial epoch 
corresponds to the phase of the maximum of the modula- 
tion. It can be proved analytically that the phases of the 
components of a modulation frequency-triplet are identical 
if the modulation is pure amplitude modulation and the ini- 
tial epoch corresponds to 7r/2 using sine terms, i.e, it is cho- 
sen t o be at the maximum phase of the modulation (|Szeidll 
120091 ). The phase coherency of the triplets and the symme- 
try of the amplitudes of the side frequencies confirm that 
the modulation of the time transformed data is basically 
amplitude modulation. 

Base d on analytic al and test results it was recently 
shown by ISzei dj (|2009h that the phase difference between 
the phases of the amplitude and phase modulations is con- 
nected to the difference between the squares of the ampli- 
tudes of the fo + fm and fo — fm components. If the f + f m 
component has larger amplitude than fo — f m has, then the 
occurrence of the largest amplitude of the pulsation (max- 
imum of the brightness maxima) precedes the occurrence 
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Figure 8. Fourier amplitudes of the first 5 harmonic components 
of the V light curves (left-hand panels) and the phase of the /o 
pulsation frequency and the epoch independent phase differences 
(right-hand panels) at 10 different phases of the Blazhko cycle. 
Note, that the scales of each amplitude and phase plots are iden- 
tical. 
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Figure 9. Folded light curve of the V observations (left-hand 
panel) and that of the time transformed data (right-hand panel). 
The times of the observations are corrected according to the phase 
variation of the fo pulsation component during the modulation 
cycle in the transformed data. The inserts magnify the plots at 
around the middle of the rising branch. The original data show 
significant spread of the observations here due to phase modula- 
tion, while in the time transformed data the phase modulation is 
eliminated as the narrowness of the rising branch indicates. 
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Figure 10. Phases of the pulsation (fc/o) and modulation side 
lobe (fc/o ± fm) frequencies for the original and the time trans- 
formed V data sets are shown in the top panels. The phase differ- 
ences between the three components of the triplets remain close 
the same in the different orders with larger deviations from these 
constant values in the k = 1, 2 and in the highest orders in the 
original V data (left-hand panel). The phases of the components 
of the triplets have nearly the same values in the different or- 
ders in the time transformed data (right-hand panel). Though 
the phases of the different frequency components are epoch de- 
pendent, the rough constancy of the phase differences between 
the triplet components at different orders holds for any epoch 
both for the original and the time transformed data. According 
to model calculations, the phase differences are close to zero (as 
observed in the case of the time transformed data) if only ampli- 
tude modulation occurs and the initial phase corresponds to the 
maximum amplitude of the modulation. Bottom panels show the 
amplitudes of the side lobe frequencies for the original and the 
time transformed data. The fc/o + fm components have systemati- 
cally larger amplitudes than the fc/o — fm components have in the 
original data as shown in the bottom left panel. Contrarily, the 
amplitudes of the side frequencies are very similar in each order 
in the time transformed data (bottom right panel), as expected 
if the modulation is pure amplitude modulation. 



of the largest delay of the maximum timings (largest pos- 
itive O — C value) and the direction of the progression in 
the maximum brightness-maximum light phase plot (right- 
hand panel in Fig. [3]) is anti-clockwise. Contrarily, if the 
larger amplitude modulation components are the fc/o — f m 
frequencies, then the phase of the largest delay of the light 
curve precedes the phase of the maximum amplitude. DM 
Cyg is an example for the former case, the amplitudes of the 
fc/o + fm components are larger than the amplitudes of the 
fc/o — fm components, the maximum of the amplitude modu- 
lation precedes the phase of the largest delay of the maxima 



and the direction of the progression in the right-hand panel 
in Fig. [3] is anti-clockwise. 

Due to the non-sinusoidal shape of the light curve and 
the residual phase modulation in the higher order pulsation 
components the phases of maxima vary slightly in the time 
transformed data as well. The difference between the phase 
of the highest light maximum and the phase of the maximum 
brightness of the largest delay during the Blazhko cycle is 
around 180° for this data set as it is shown in Fig. 1111 i.e. the 
largest delay of the light curve occurres when the amplitude 
of the pulsation is the smallest. Fig.[TT]shows the same plots 
as Fig. [3] but for the time transformed data. The right-hand 
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Figure 11. The same figure as Fig. [3] but for the time transformed data. The phase difference between the maximum brightness and 
maximum phase timings for the time transformed data is around 180°. The maximum brightness - maximum phase plot has been 
substantially compressed, the direction of the progression is denoted by an arrow. Dashed lines connect the mean values in 10 phase bins 
of the modulation. 



panel shows already hardly any loop structure wider than 
the scatter of the data. 



3.2 Changes in the global physical parameters 
during the Blazhko cycle 

In ISodor. Jurcsik fc Szeidll (|2009h we demonstrated that 
from good quality multicolour photometric data the physi- 
cal parameters of RRab stars can be determined with similar 
accuracy as with direct Baade-Wesselink analysis, using an 
inverse photometric method (IPM) . We applied this method 
to the BVIc light curves of MW Lyr at different phases of 
its Blazhko modulation succesfully, about 1 — 2% changes 
in the mean global physical para meters (L, T e f f , R) were 
detected during the Blazhko cycle (I Jurcsik et alj|2009l ). 

Before analysing the light curves of DM Cyg at different 
phases of its modulation, first we have to determine those 
parameters that definitely do not vary during the Blazhko 
cycle. These are the metallicity, the mass and the distance 
of the star and the dereddened standard magnitudes of the 
comparison stars i.e., the zero points of the magnitude scales 
used. 

The [Fe/H] of DM Cyg w as det ermined spectr oscop- 
ically by ISuntzeff et al.l (|l994f ) and iLavdenl (|l994h . The 
Fourier paramete rs of the mean V l i ght cu rve give [Fe/H] = 
—0.01 using the lJurcsik fc Kovacsi (|l996l . Eq. 3) formula. 
As the spectroscopic observations correspond to —0.16 
and 0.07 [Fe/H] values on the metallicity scale of the 
lJurcsik fc Kovacsi (|l996h formula, atmosphere models with 



[Fe/H ]= 0.0 and —0.1 compositions (|Castelli fc Kuruczl 
120031 ) are used in the light curve modelling of DM Cyg. 

The IPM finds mass values that are too large on evolu- 
tionary grounds in some cases if the mass is also allowed to 
vary by the fitting process. The fitting accuracy of the dif- 
ferent mass solutions differs, however, only marginally. The 
mean B,V,Ic light curves of DM Cyg are fitted with 5.1, 
4.7 and 4.0 mmag residual scatters when the mass value is 
also allowed to vary. In this case the best solution is found at 
m = 0.84SHs un . The r.m.s. of the fixed M = 0.55501^ mass 
solution B, V, Ic light curve fits are 5.2, 4.6 and 4.3 mmag, 
respectively. The average fitting accuracy in the three bands 
decreases only 0.1 mmag, from 4.7 to 4.6 mmag, when the 
mass is also allowed to be fitted. That is, the method is 
rather insensitive to the value of the mass. Therefore, in or- 
der to obtain reliable solutions, the IPM has been run with 
0.50, 0.55 and 0.60 9Jlsun fixed mass values for the mean 
light curves of DM Cyg. 

The global mean absolute physical parameters of DM 
Cyg averaged over both the pulsation and the Blazhko cycles 
and its distance are derived from the mean light curves using 
the IP method. The results are summarized in Table [3] as- 
suming different possible mass and [Fe/H] values. Note that 
the zero points of the colour and magnitude scales defined 
by the dereddened magnitudes of the comparison stars, are 
also allowed to vary i.e, they are also fitted at this step. To 
derive the distance, the apparent dereddened mean V mag- 
nitude of DM Cyg has to be known as the IPM gives the 
absolute mean V brightness of the star. To determine the 
apparent dereddened mean magnitude of DM Cyg, the stan- 
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Table 4. Mean physical parameters of DM Cyg derived from its mean light curves using the IP method. 
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Figure 12. Variations of the observed mean (left-hand panels) and derived parameters of DM Cyg during the Blazhko cycle. Circles and 
crosses denote the results when the IP method use Liu's template radial velocity curve, and when the V ra d curve is calculated from the 
Iq light curve, respectively. No significant, systematic difference between the results can be seen depending on the choice of the initial 
radial velocity templates. The results of the IP method for fixed 0.50, 0.55 and 0.60 OT/3Jt Slln mass values show the same variations in 
each parameters. Only the absolute mean values of the parameters are different, according to the data given in Table [4] These results 
prove that the detected changes of the physical parameters during the Blazhko cycle of DM Cyg are not sensitive to the settings of the 
IP method, and are independent of the uncertainties of the values of the mean global parameters. See further details in the text. 



dard V magnitude of the comparison stars and the interstel- 
lar absorption Ay should have to be known. Unfortunately, 
no standard magnitudes of the comparison stars have been 
published. We estimate the average of the dereddened V 
brightnesses of the two comparison stars to be 11.333 mag 
trans formed from their Tycho Bt, and Vt magnitudes <|ESAl 
1 19971 ) and using A v = 3.UE(B -V) = 0.558 interste llar ab- 
sorption value according to the ISchlegel et al.1 (|l998l ) maps. 
It has to be emphasized, however, that the uncertainties of 
the standard V magnitudes of the comparison stars and the 



interstellar absorption value affect only the distance esti- 
mate. Neither any other parameter nor their detected varia- 
tions depend on the choice of the dereddened V magnitudes 
of the comparison stars. 

The uncertainties of the estimates of the absolute phys- 
ical parameters listed in Table [4] correspond to the standard 
deviations of the results of running the IP code with 16 
different settings. These settings are the combinations of: 

i) B, V and Ic, or V, B — V and V — Ic data are used; 

ii) the V ra d curve is defined by Liu's template or it is cal- 
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culated from the Ic light curve; Hi) 2 different weights of 
the initial V ra d curve are used; iv) 2 different values of the 
AA(y ra d) I AVAmp ratio valid for Blazhko variables are ap- 
plied (see further d etails in ISodor. Jurcsik fc Szeidll |2009| ; 
I Jurcsik et aTll2009l ). However, as it can be seen from the 
data listed in Table [3] the true uncertainties of the mean 
global parameters are larger than the estimated errors of 
the solutions obtained for fixed mass and metallicity values. 
The true possible parameter space comprises the full range 
of the solutions for the possible mass and metallicity ranges. 
Nevertheless, we are focusing on the variation in the mean 
physical parameters and not on their absolute values, and it 
is shown in the rest of this section that the uncertainties of 
the mean global physical parameters have no effect on their 
variation during the Blazhko cycle. 

Though the amplitude of the modulation of DM Cyg is 
only about one tenth of that of MW Lyr, we have applied 
the IPM to the light curves in 10 different phases of the mod- 
ulation in order to decide whether there are any detectable 
changes in the mean global parameters of DM Cyg during 
its Blazhko period. The distance and the mass are fixed to 
their possible values given in Table [4] assuming [Fe/H]= 0.0 
metallicity. The zero points of the magnitude scales are also 
fixed correspondingly to the solution obtained for the mean 
light curves. 

Fig. [12] shows the results for the solutions using 9JI = 
0.509Jt Su n, d = 1179 pc; DJl = 0.55OT Stl n, d = 1215 pc and 
QTt = 0.609Jtsun, d = 1233 pc mass and distance values. In 
the left-hand panels, the variations in the observed mean 
quantities are plotted: total pulsation amplitude in V band, 
the variation of the pulsation period derived from the phase 
variation of the /o pulsation frequency and different pulsa- 
tion averages of the magnitudes and colours. The right-hand 
panels show the quantities derived from the IP method for 
three possible mass/distance combinations: the mean val- 
ues of the radius, the surface gravity, the absolute visual 
brightness averaged by magnitude and intensity units, the 
luminosity and the effective temperature. The different com- 
binations of the mass and distance give very similar results, 
the amplitudes and phases of the detected changes in the 
mean physical parameters during the Blazhko cycle hardly 
change, only their averages over the Blazhko cycle are dif- 
ferent, corresponding to their respective values given in Ta- 
ble[4] If atmosphere models with [Fe/H] =— 0.1 are used, the 
results on the variations of the physical parameters during 
the Blazhko cycle are not changed. 

In order to get a real estimate of the uncertainties of 
the derived quantities in different phases of the modulation, 
the IP code has been run with 16 different settings for each 
data set again. Note that the zero points of the colour and 
magnitude scales defined by the dereddened magnitudes of 
the comparison stars, are also fixed now. These magnitudes 
were, however, allowed to vary i.e., they were also fitted 
when the mean global parameters were determined. This 
explains why the scatter of the results for the 16 different 
setting of the IP method is much smaller in the different 
phases of the modulation than the scatter of the results for 
the mean light curves. 

The different intensity and magnitude averages of the 
V light curve and the colour curves show hardly any vari- 
ations, the amplitudes of their changes is only 0.002-0.004 
mag. The V, and B — V averages show some small system- 



atic variations of similar character like these averages in MW 
Lyr, but the V — I averages show only scatter. In spite of 
that the observed mean magnitudes and colours vary only 
slightly during the Blazhko cycle in DM Cyg, by the aid 
of the IP method 0.3% and 7 K systematic changes in the 
mean luminosity and temperature of the star could be de- 
tected during its Blazhko cycle. The phase relation of these 
variations are the opposite, similarly as it was detected in 
MW Lyr. Both stars are the most luminous and the coolest 
at around the largest amplitude phase of the modulation. 

The IP method calculates the mean physical parame- 
ters as their mathematical averages over the pulsation cycle. 
F($) =< F(<p) >, where $ and tp denote modulation and 
pulsation phases, respectively, and <> denotes averaging by 
ip for the whole pulsation cycle. For the luminosity, it means 
that the values plotted in Fig. 1 121 correspond to: 

L($) = < Hip) > = < 4naR(ip) 2 T eff {ip) A >, 
as the IP method satisfies the Stephan-Boltzmann law in 
eac h phase of the p ulsatio n (ip). 

ICarnev et all (|l992l ) showed from direct Baade- 
Wesselink analysis results that, even when the equilibrium 
luminosity and radius of pulsating variables equal with the 
arithmetic means of their variations over the pulsation cy- 
cle, the equilibrium temperature T eq is not the same as the 
mean temperature < T e // > of the star. The equilibrium 
temperature is defined then as 

T eq = (L eq /4naR 2 eq )^ 4 = (L($)/4^aJ?($) 2 ) 1 / 4 
= (< L{<p) > /47ia < R{ip) > 2 ) 1/4 . 
Consequently, the relation between the arithmetic means of 
the physical parameters shown in Fig. [12] deviates to some 
extent from the theoretical Stefan-Boltzmann law. 

The 0.13% variation in the mean radius is in good agree- 
mant with the 0.17% changes of the pulsation period (the 
pulsation equation requires AP/P ~ 3/2AR/R). The pul- 
sation period changes are determined from the phase differ- 
ences of the /o pulsation frequency in different phases of the 
modulation, but when deriving the radius variation the IP 
method does not utilize this information in any way. Con- 
sequently, the determined period and radius variations are 
completely independent quantities. The fact, that they show 
such a good agreement is a great support both to the IP 
method, and to our interpretation of the phase modulation 
as variations in the pulsation period. 

Comparing the res ults obtained for DM Cyg with that 
of MW Lyr (Fig. 14 in I Jurcsik et alj|2009h the similar be- 
haviour of the two stars are conspicuous. The only difference 
is in the detected amplitudes of the variations of the differ- 
ent parameters, each being about one tenth of their detected 
amplitudes in MW Lyr, in accordance with the different am- 
plitudes of the modulation of the two stars. 



4 THE PHOTOGRAPHIC AND 
PHOTOELECTRIC DATA 

The p hotographic light curve of DM Cyg was shown in lHurtal 
(2009). Although the Konkoly photographic data are very 
sparse, the observations covered 24 years with a 13-year gap, 
we tried to analyse this data set also searching for any sign 
of light curve modulation. No indication of light curve mod- 
ulation has been found in any of the two parts of the pho- 
tographic data, but there is some hint that modulation fre- 
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Figure 15. The O — C diagram of DM Cyg. Literature data are 
from the GEOS database (dots), the most deviant data are omit- 
ted. The Konkoly photographic, photoelectric and CCD normal 
maximum timings listed in Tabic [2] are plotted by squares, the 
time intervals of the Konkoly observations are indicated by hori- 
zontal lines in the top of the figure. The O — C indicates steady 
period increase during the ~ hundred years of the observations. 



Figure 13. Amplitude spectrum of the residual photographic 
light curve (after the removal of the pulsation signal) in the vicin- 
ity of the /o,2/o and 3/o pulsation frequencies. The residuals 
show signals at around /o + f m , 2/o + f m and 3/o + fm frequen- 
cies with 0.02 — 0.03 mag amplitude. The fo and f m frequencies 
detected in the photographic data are 2.3817564 and 0.0940, re- 
spectively. Bottom panel shows the spectral window of the pho- 
tographic observations. 
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Figure 14. Maximum brightness and phase of the B pg = 
11.5 mag brightness values on the rising branch (0.2 mag below 
maximum brightness) for the photographic observations phased 
with the supposed modulation period (10.636 days) are shown in 
the left-hand and right-hand panels, respectively. Data from the 
two parts of the photographic observations are shown with dif- 
ferent symbols. According to these plots, there is some indication 
that the maximum brightness and the phase of the rising branch 
of DM Cyg varied with small amplitude during the time of the 
photographic observations as well. 



quency components appear at fc/o + fm, (k = 1, 2, 3) posi- 
tions in the residual spectrum of the complete photographic 
data set (see Fig. I13[l . The highest modulation peaks are not 
exactly at the same separations for fo + f m , 2/o + f m and 
3/o + f m - The best solution which gives the smallest resid- 
ual for the photographic data can be gained with 2.3817564 
c/d (0.41985822 d) and 0.0940 c/d (10.636 d) pulsation and 
modulation frequencies, respectively. 

Fig. [14] plots the maximum brightness and the phase 



phased with the supposed 10.636 d modulation period. The 
phase of a given magnitude on the rising branch can be de- 
termined with higher accuracy than the phase of the max- 
imum, therefore we use this quantity to measure the am- 
plitude of the phase modulation component in the photo- 
graphic data. Data belonging to the two parts of the obser- 
vations are denoted by different symbols in Fig. Q3] Accord- 
ing to these plots the amplitudes of the amplitude and phase 
modulations were not larger at the time of the photographic 
observations than today. 

The photographic data were also analysed taking into 
account the period decrease of the pulsation that took place 
during the time interval of the observations, but the results 
were only marginally different from that obtained from the 
original photographic data. 

The photoelectric observations comprise data only from 
four nights which does not allow to check the light curve 
modulation. The photoelectric data were therefore used only 
to determine one normal maximum timing listed in Table [2] 

The normal maximum timings of the Konkoly obser- 
vations given in Table [2] together with literature data col- 
lected in the GEO S database are plotted in Fig. 1151 Based 
on the GEOS data lLe Borgne et all (|2007T ) determined p = 
0.091d Myr -1 period increase rate for DM Cyg. The addi- 
tion of the Konkoly data does not modify this result. 

From the analyses of the CCD and photographic data of 
DM Cyg we conclude that the period changes of the pulsa- 
tion and modulation have the opposite direction, when the 
pulsation period of DM Cyg was shorter by 0.000005 d then 
the modulation period was 0.066 d longer than today. These 
values correspond to ^ = -13200 or ^pf/^ = -523 pe- 
riod change rates. However, these data have to be taken with 
caution partially because the uncertainties of the results de- 
rived from the photographic data and also because pulsation 
and modulation period changes of those Blazhko variables 
where the periods could be determined for several epochs 
show that there is no strict relation between the two periods, 
for some time intervals even the sign of the period change 
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ratio can change (e.g., in RV UMa lHurta et al" . 2008]). There- 
fore, the period change ratios defined by only two-epoch data 
should be missleading in some cases. We also have to note, 
however, that the complex period change behaviour of the 
modulation is usually connected to complex changes of the 
pulsation period and of the pr operties of the modulation, too 
(e.g., in RR Gem and XZ Cve lSodor, Szeidl fe Jurcsikll2007l ; 
lLaCluvze et alj|2004h . Probably, the steady period change 
of the pulsation of DM Cyg is accompanied with steady pe- 
riod change of its modulation and with the stability of its 
modulation properties. 



5 SUMMARY 

In the present paper the light curve modulation of DM Cyg, 
a fundamental mode RR Lyrae variable has been investi- 
gated. Most of the results are very similar to those found 
in the analysis of MW Lyrae, a large modulation amplitude 
Blazhko variable, but on a significantly smaller scale. The 
0.07 mag amplitude of maximum brightness variation of DM 
Cyg is only 10% of the modulation amplitude of MW Lyr, 
accordingly, all the changes detected in DM Cyg during its 
Blazhko period are about 10% of the detected changes in 
MW Lyr. 

The phase and amplitude relations of the amplitude 
and phase modulation components are similar in DM Cyg 
and MW Lyr. The variations of the mean magnitudes and 
colours in DM Cyg where they can be detected are also sim- 
ilar to their counterparts in the MW Lyr data. Therefore, it 
is not surprising that the derived changes in the mean global 
physical parameters of DM Cyg resemble the changes found 
in the mean parameters of MW Lyr only on a much reduced 
scale. 

The only notable difference between the character of 
the modulations of the two stars is found when the changes 
of the Fourier parameters of the light curves during the 
Blazhko cycle is investigated. While in the case of MW Lyr 
the phases of 2/o, and 3/o change in line with the phase vari- 
ation of /o (i.e, there is no significant changes detected in the 
f(fki) phase differences), when analysing the light curves of 
DM Cyg 90° phase differences between the phase variation 
of the /o pulsation frequency and the phase variations of 
the higher order pulsation components are found. In spite of 
this, the data transformation that corrects the times of the 
observations taking into account the phase variation of the 
/o pulsation frequency which completely separated the am- 
plitude and phase modulation components of the light curve 
modulation of MW Lyr seems to work also on the light curve 
of DM Cyg. 

The modulation behaviour of other Blazhko variables 
can differ, however, more significantly from the modulations 
of DM Cyg and MW Lyr. There are Blazhko stars, where the 
phase relation of the amplitude and phase modulation com- 
ponents contrast with the phase relations of DM Cyg and 
MW Lyr. Also, the detected changes in the mean colour s 
may vary differently e.g, in SS Cnc (| Jurcsik et al.l l2006h . 
Therefore, we cannot draw major conclusions from the sim- 
ilarity of the results obtained for DM Cyg and MW Lyr be- 
fore analysing other Blazhko variables of different character 
on a similar way. 
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ABSTRACT 

DM Cyg, a fundamental mode RRab star was observed in the 2007 and 2008 
seasons in the frame of the Konkoly Blazhko Survey. Very small amplitude light curve 
modulation was detected with 10.57 d modulation period. The maximum brightness 
and phase variations do not exceed 0.07 mag and 7 min, respectively. In spite of the 
very small amplitude of the modulation, beside the frequency triplets characterizing 
the Fourier spectrum of the light curve two quintuplet components were also identified. 
The accuracy and the good phase coverage of our observations made it possible to 
analyse the light curves at different phases of the modulation separately. Utilizing the 
IP method (Sodor. Jurcsik and Szeidl, 2009) we could detect very small systematic 
changes in the global mean physical parameters of DM Cyg during its Blazhko cycle. 
The detected changes are similar to what we have already found for a large modulation 
amplitude Blazhko variable MW Lyrae. The amplitudes of the detected changes in the 
physical parameters of DM Cyg are only about 10% of that what have been found in 
MW Lyr. This is in accordance with its small modulation amplitude being about one 
tenth of the modulation amplitude of MW Lyr. 

The pulsation period of DM Cyg has been increasing by a rate of fj = 0.091dMyr _1 
during the hundred-year time base of the observations. Konkoly archive photographic 
observations indicate that when the pulsation period of the variable was shorter by 
Appuls = 5 • 10~ 6 d the modulation period was longer by Ap nlo d = 0.066 d than today. 

Key words: stars: variables: other - stars: horizontal branch - stars: individual: 
DM Cyg - techniques: photometric - methods: data analysis - 



1 INTRODUCTION 

Utilizing our full access to an automatic 60 cm telescope we 
have obtained extended multicolour observations of many 
fundamental mode RR Lyrae variables showing light curve 
modulation (the Blazhko effect) during the past five years. 
Detailed analyses of some of our targets were already pub- 
lished in I Jurcsik et~all (|2005l , 120061 , 120081 . 120091 ). Our obser- 
vations are the first multicolour photometric data which are 
accurate and dense enough to allow not only to determine 



* Based on observations collected with the automatic 60 cm tele- 
scope of Konkoly Observatory, Budapest, Svabhegy 



the frequencies appearing in the Fourier spectra of the light 
curves, but also to find changes in the global mean physical 
parameters of the stars (L,T e ff,R) in different phases of 
the modulation. In order to extract changes in the phys- 
ical parameters exclusively from multicolour photometric 
data we devel oped an inverse photometric B aade-Wesselink 
method (IPM: ISodor. Jurcsik fc Szeidlll2009l ). The IPM was 
successfully applied to determine changes in the mean global 
physical parameters of MW Lyr during its Blazhko cycle in 
lJurcsik et all (|2009h . 

DM Cyg, a short period (P = 0.42 d) RRab star was 
selected to be observed in the frame of the Konkoly Blazhko 
Survey because it had been announced to show phase modu- 
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Table 1. Konkoly photometric observations of DM Cyg. 



HJD 


mag* 


detector 


filter 


2400000+ 








54265.45479 


-0.249 


CCD 


V 


54265.46045 


-0.404 


CCD 


V 


54265.46611 


-0.501 


CCD 


V 


* Relative magnitudes for the CCD and photoelectric data. 



Table 2. Normal maximum timings derived from the Konkoly 
photographic, photoelectric and CCD observations. 



.ID 



Normal maximum timings obs. 



2427667 - 2429226 
2434186 - 2436373 
2443747 - 2444133 
2454265 - 2454407 
2454597 - 2454711 



2429078.482 

2435342.347 

2443777.333 

2454334.3467 

2454661.4210 



Pg 

Pg 

pe 
CCD 
CCD 



lation with a period of 26 d by Lvsova fc Firmanvukl (|l980t ). 
but neither the NSVS data jWozniak et all bOO^ Tnor the 
complete data set of its maximum timings show phase mod- 
ulation with this period (|S6dor fc Jurcsikll2005l ). 

Though DM Cyg is a relatively bright (V = 10-11 mag) 
RR Lyrae star, complete, accurate, multicolour light curve 
of its pulsation has never been published. Its period change 
was, however, regularly monitored by different groups of ob- 
servers. The GEOS databastf] lists 260 maximum times of 
DM Cyg between 1900 and 2008. 

In the present paper we publish our extended CCD ob- 
servations of DM Cyg and the results of the analysis of its 
light curve modulation. Archive photographic and photo- 
electric Konkoly data are also processed. 



2 DATA 

CCD observations were obtained with the automated 60 
cm telescope of the Konkoly Observatory, Svabhegy, Bu- 
dapest equipped with a Wright Instruments 750 x 1100 
CCD camera and BVIc filters. Measurements were taken 
on 81 nights between July 2007 and Sept 2008. About 
3100 data points in each band were gathered. Exposition 
times were 200, 60 and 40 sec or a bit longer depend- 
ing on the sky transparency in the BVIc bands, respec- 
tively. Data reduction was performed using standard IRAF0 
packages. Aperture photometry of DM Cyg (21:21:11.548 
+32:11:28.71) and several neighbouring stars were carried 
out in order to check the stability of the photometry and 
the constancy of the comparison stars. The relative mag- 
nitudes of DM Cyg measured to the mean magnitudes of 
C1=GSC2.2 N0330220980 (21:21:30.731 +32:13.05.78), and 



1 http : //dbrr . ast . obs-mip . f r /maxRR . html 

2 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Figure 1. Delta V magnitudes versus Julian Date, data phased 
with the 0.419863 d pulsation and the 10.57 d modulation periods 
are shown in panels A, B and C, respectively. 



C2=GSC2.2 N03302207371 (21:20:59.483 +32:13.01.63) are 
used in the analysis. The Tycho B and V magnitudes of the 
comparison stars CI and C2 are B T = 12.705, V T = 11.799, 
and B T = 13.081, V T = 12.150, respectively (|Hog et all 
|2000| ). Second order extinction correction of the data were 
applied in the B band. Relative magnitudes are transformed 
to standard BVIc magnit udes. More details o n the reduc- 
tion procedure are given in I Jurcsik et al.l (|2008l ). 

Archive photoelectric and photographic data obtained 
with the 60 cm telescope in 1978 and with a 16 inch astro- 
graph between 1934 and 1958 are also utilized. Photoelectric 
observations were obtained on 4 nights, the BV magnitudes 
were measured relative to GSC2.2 N03302207371. The pho- 
tographic measurements comprise data from 40 nights, the 
plates were evaluated using GSC2.2 B magnitudes of the 
surrounding stars. 

Table [1] shows a sample of the Konkoly photometric 
data of DM Cyg. All the photometric observations are avail- 
able in the online version of the journal as Supplementary 
Material. 

Seasonal normal maximum timings from each data set 
have been determined for five epochs, these data are listed 
in Table [2] 



3 RESULTS 

The accurate CCD observations reveal that the light curve of 
DM Cyg is not stable. It shows small amplitude modulation, 
the variations in maximum brightness and phase are only 
about 0.07 mag and 0.005 d (7 min, 0.012 pulsation phase), 
respectively. 

The V light curve of DM Cyg and the results obtained 
folding the data with the pulsation and modulation periods 
are shown in the three panels of Fig. [T] The phased resid- 
ual light curve after the removal of the mean pulsation light 
variation is given in Fig. [2] This figure shows the differences 
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Figure 3. Maximum V brightness and maximum phase values phased according to the Blazhko period are plotted in the top-left and 
bottom-left panels, respectively. Their typical errors are 0.005 mag and 0.002-0.004 phase. The right-hand panel shows the maximum 
brightness vs maximum phase data. The progression of data during the Blazhko cycle in this diagram is anti-clockwise. Dashed lines 
connect the mean values in 10 phase bins of the modulation. 
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Figure 2. Residual V light curve of DM Cyg after removing the 
pulsation components from the data according to the light curve 
solution given in Table \3\ Residuals of the smallest (stars) and 
the largest (diamonds) amplitude phases are shown with large 
symbols in the left-hand and right-hand panels, respectively. Note 
that the residuals have larger amplitude at around the middle of 
the rising branch of the pulsation (phase = —0.1) than around 
pulsation maxima (phase = 0) due to phase modulation. 

of the amplitude of the modulation at different phases of the 
pulsation. The well observable amplitude of the modulation 
is concentrated to a narrow phase range of the pulsation 
around minimum, rising branch and maximum. A similar 
behaviour has been already found in RR Gem and SS Cnc, 
in two other small modula tion amplitude Blazhko variables 
jjurcsik et al.ll2005l , [200(3 ). The amplitude of the the resid- 



ual variation is the largest at around the middle of the rising 
branch of the pulsation light curve (phase = —0.1) indicat- 
ing phase modulation of the rising branch as well. There is a 
difference between the phase of the amplitude modulation of 
the maximum of the light curve and the phase of the phase 
modulation of the same feature. In particular the maximum 
positive displacement of the timings of the maximum bright- 
ness precedes the occurrence of the brightest maximum by 
about 0.25 Blazhko phase as it is shown Fig. [3] 

The elements of the pulsation and the modulation are: 

T maxpuls = 2 454 312.514 [HJD] + 0.419863 ■ E puls , (1) 

and 

TmaxBi = 2 454 312.514 [HJD] + 10.57 ■ E m . (2) 

-Epuls and Ebi denote the epoch number of the pulsation and 
modulation cycles, respectively. 

The pulsation and modulation periods correspond to 
the mean values of the frequencies of the light curve so- 
lutions of the B, V and Ic data using 'locked' frequency 
solutions of the kfo pulsation and fc/o ± / m , and / m modu- 
lation frequency components (i.e., the modulation side lobe 
components are at the positions of the linear combination 
frequencies) . 

Data analysis was performed u sing the differ ent appli- 
cations of the MUFRAN package l |Kollath|[i99ch . a linear 
combination fitting program developed by A. Sodor, and 
the linear and nonlinear curve fitting abilities of gnuplo10. 

3 http://www.gnuplot.info/ 
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Figure 4. The left-hand panels show the spectral window, the amplitude spectrum and the residual spectrum after the removal of the 
pulsation components for the CCD V data. The middle and right-hand panels show residual spectra in the vicinity of the 10/o,H/o, 
and 12/o (fc = 10, 11, 12) and of the 25/o, 26/o, and 27/o (k = 25, 26, 27) frequencies, respectively. The spectrum has been prewhitened 
for the pulsation components up to the 24th order, and in the right hand panels for the detected modulation components, too. These 
spectra indicate that modulation frequency components are present only in the first 10-11 orders, while the pulsation can be accurately 
described if harmonic frequency components are taken into account up to as high as the 27th order. For clarity, vertical grids denote the 
positions of the fc/o, fc/o + fm and kfo — f m frequencies. 
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Figure 5. Fourier spectra of the V and B residual data in 
the vicinity of /o and 2/rj. The pulsation and the fc/o ± f m 
modulation frequencies have been removed. The /o — 2/ m and 
2/o — 2fm quintuplet components appear in these residual spec- 
tra with ~ 0.001 mag amplitude. 
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3.1 The light curve solution 

The Fourier amplitudes and phases of the pulsation and 
modulation frequency components identified in the spectra 
of the BVIc light curves of DM Cyg are summarized in 
Table [3] The Fourier decompositions use sin terms and the 
initial epoch corresponds to one of the brightest maxima of 
the modulated light curve, TO = 2454312.514. The errors of 
the amplitudes are ~ 0.0002 mag, the errors of the phases for 
the time transformed V data (see the details later) are given 
in the last column. The pulsation components are detectable 
up to the 27th order, while the modulation side frequencies 
(fc/o±/m) diminish at around the 11th order as documented 
in Fig. [4] The modulation frequency f m is present without 
question in each of the B,V,I spectra with similar ampli- 



Figure 6. Amplitude decrease of the pulsation and modulation 
side lobe frequencies at different orders shown as amplitudes nor- 
malised to the amplitudes of the first order (k=l) components. 
The decrease of the amplitudes of the harmonic components of 
the pulsation is exponential like, while the decrease of the ampli- 
tudes of the side lobe frequencies is more linear, with irregular 
character in the low orders. 



tudes as the 4th order negative and the 9th order positive 
modulation side lobe components have. Quintuplet frequen- 
cies are detected at fo — 2f m and 2/o — 2/ m as shown in 
Fig. El 

The 51 frequencies listed in Table[3]fit the data with 6-8 
mmag residual scatter in each band, which is about the level 
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Table 3. Fourier amplitudes and phases of the pulsation and modulation frequencies in DM Cyg. Initial epoch corresponds to maximum 
of both the pulsation and the modulation light variations as given in Eqs. 1 and 2. 



frequency B V I c V t * t a v{V tt ) 
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0.0015 
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* Fourier parameters of the time transformed V data 



of the observational noise. No further periodic or stochastic 
variation in the residuals is detected. 

The decrease of the amplitudes of the pulsation compo- 
nents with increasing order is smooth and exponential-like. 
The amplitude decrease of the modulation components is 
however, different. The amplitudes of the low order mod- 



ulation components behave irregularly, while they show a 
linear decrease at higher orders (see Fig. IfoJ). Similar be- 
haviour of the amplitude decrease of the detected frequen- 
cies were demonstrated for RR Gem, SS Cnc a nd SS For 
(IJurcsik et al. I l2005l . 12006] ; iKolenberg et al.ll200Sl ) . 

Fig.[7]shows the V light curves belonging to 10 different 
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Figure 7. AV light curves in 10 phase bins of the modulation. All the observations are also shown with gray colour in each plot. Empty 
circles denote artificial data which are used to stabilize the Fourier fits if there is some gap in the observations. 



phase bins of the modulation. The first five Fourier param- 
eters of the 15th order harmonic fits to these individual V 
light curves are plotted in Fig. [8] Comparing the variations 
of the Fourier parameters of the l ight curves of DM Cyg 
with that of MW Lyr (Fig. 12 in I Jurcsik et al.l l2008h the 
following differences are conspicuous. The amplitudes of the 
detected changes in the Fourier parameters are about one 
order of magnitude smaller than the amplitudes of the light 
curve variations in MW Lyr. The amplitudes of the varia- 
tion of the amplitudes of the fo, ...5/o components decrease 
more drastically with increasing order in MW Lyr than in 
DM Cyg. This is most probably connected to that the pul- 
sation components are detected only up to the 12th order in 
MW Lyr, while they can be observed up to the 27th order in 
DM Cyg. Nevertheless the phase relations between the am- 
plitude modulation and the variation in the phase of the fo 
pulsation component are the same for the two stars, the de- 
tected changes in the epoch independent phase differences 
[v(/fei)] are, however, significantly different. While in DM 
Cyg the tp(fki) phase differences show sinusoidal variations 
with 0.25 (it/2) phase shift and with larger amplitude than 
<p(fo), in MW Lyr the amplitudes of the variations of the 
P(fki) components are smaller than the amplitude of <p(fo) 
and they do not show strictly regular behaviour with the 
modulation period. 

In lJurcsik et al l i|200Sl ) we have introduced a new 
method of analysing Blazhko variables' light curves. It was 
shown that if the times of the observations are corrected 
according to an appropriate time transformation defined by 
the variation of the phase of fo during the Blazhko cycle, 
then the modulation can be separated into phase and am- 
plitude modulation components. The time transformed data 
showed pure amplitude modulation for MW Lyr. Though the 
modulation amplitude of DM Cyg is very small, we have 
checked how a similar time transformation influences its 
modulation properties. Transforming the times of the ob- 



servations according to a continuous harmonic function of 
the tp(fo) data, the V observations folded with the pulsation 
period is compared to the folded light curve of the original 
data of DM Cyg in Fig. Enlarged plots of the middle of 
the rising branch are inserted in the figures. Though the dif- 
ference between the original and the time transformed data 
is small, the reduction of the phase modulation component 
in the time transformed data is indicated by the narrowness 
of the rising branch in the right-hand panel of Fig. [5] 

In the last three columns of Table [3] the Fourier ampli- 
tudes and phases and the errors of the phases are listed for 
the time transformed V data. The amplitudes of the mod- 
ulation side frequency components kfo + f m a,nd kfo — f m 
are about the same and the phases of the triplets in the dif- 
ferent orders show definite phase coherency according to the 
Fourier solution of the time transformed data (see Fig. 1 10p . 
As the phases of the components of the triplets are epoch de- 
pendent, this phase coherency holds only if the initial epoch 
corresponds to the phase of the maximum of the modula- 
tion. It can be proved analytically that the phases of the 
components of a modulation frequency-triplet are identical 
if the modulation is pure amplitude modulation and the ini- 
tial epoch corresponds to 7r/2 using sine terms, i.e, it is cho- 
sen t o be at the maximum phase of the modulation (|Szeidll 
120091 ). The phase coherency of the triplets and the symme- 
try of the amplitudes of the side frequencies confirm that 
the modulation of the time transformed data is basically 
amplitude modulation. 

Base d on analytic al and test results it was recently 
shown by ISzei dj (|2009h that the phase difference between 
the phases of the amplitude and phase modulations is con- 
nected to the difference between the squares of the ampli- 
tudes of the fo + fm and fo — fm components. If the f + f m 
component has larger amplitude than fo — f m has, then the 
occurrence of the largest amplitude of the pulsation (max- 
imum of the brightness maxima) precedes the occurrence 
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Figure 8. Fourier amplitudes of the first 5 harmonic components 
of the V light curves (left-hand panels) and the phase of the /o 
pulsation frequency and the epoch independent phase differences 
(right-hand panels) at 10 different phases of the Blazhko cycle. 
Note, that the scales of each amplitude and phase plots are iden- 
tical. 
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Figure 9. Folded light curve of the V observations (left-hand 
panel) and that of the time transformed data (right-hand panel). 
The times of the observations are corrected according to the phase 
variation of the fo pulsation component during the modulation 
cycle in the transformed data. The inserts magnify the plots at 
around the middle of the rising branch. The original data show 
significant spread of the observations here due to phase modula- 
tion, while in the time transformed data the phase modulation is 
eliminated as the narrowness of the rising branch indicates. 
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Figure 10. Phases of the pulsation (fc/o) and modulation side 
lobe (fc/o ± fm) frequencies for the original and the time trans- 
formed V data sets are shown in the top panels. The phase differ- 
ences between the three components of the triplets remain close 
the same in the different orders with larger deviations from these 
constant values in the k = 1, 2 and in the highest orders in the 
original V data (left-hand panel). The phases of the components 
of the triplets have nearly the same values in the different or- 
ders in the time transformed data (right-hand panel). Though 
the phases of the different frequency components are epoch de- 
pendent, the rough constancy of the phase differences between 
the triplet components at different orders holds for any epoch 
both for the original and the time transformed data. According 
to model calculations, the phase differences are close to zero (as 
observed in the case of the time transformed data) if only ampli- 
tude modulation occurs and the initial phase corresponds to the 
maximum amplitude of the modulation. Bottom panels show the 
amplitudes of the side lobe frequencies for the original and the 
time transformed data. The fc/o + fm components have systemati- 
cally larger amplitudes than the fc/o — fm components have in the 
original data as shown in the bottom left panel. Contrarily, the 
amplitudes of the side frequencies are very similar in each order 
in the time transformed data (bottom right panel), as expected 
if the modulation is pure amplitude modulation. 



of the largest delay of the maximum timings (largest pos- 
itive O — C value) and the direction of the progression in 
the maximum brightness-maximum light phase plot (right- 
hand panel in Fig. [3]) is anti-clockwise. Contrarily, if the 
larger amplitude modulation components are the fc/o — f m 
frequencies, then the phase of the largest delay of the light 
curve precedes the phase of the maximum amplitude. DM 
Cyg is an example for the former case, the amplitudes of the 
fc/o + fm components are larger than the amplitudes of the 
fc/o — fm components, the maximum of the amplitude modu- 
lation precedes the phase of the largest delay of the maxima 



and the direction of the progression in the right-hand panel 
in Fig. [3] is anti-clockwise. 

Due to the non-sinusoidal shape of the light curve and 
the residual phase modulation in the higher order pulsation 
components the phases of maxima vary slightly in the time 
transformed data as well. The difference between the phase 
of the highest light maximum and the phase of the maximum 
brightness of the largest delay during the Blazhko cycle is 
around 180° for this data set as it is shown in Fig. 1111 i.e. the 
largest delay of the light curve occurres when the amplitude 
of the pulsation is the smallest. Fig.[TT]shows the same plots 
as Fig. [3] but for the time transformed data. The right-hand 
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Figure 11. The same figure as Fig. [3] but for the time transformed data. The phase difference between the maximum brightness and 
maximum phase timings for the time transformed data is around 180°. The maximum brightness - maximum phase plot has been 
substantially compressed, the direction of the progression is denoted by an arrow. Dashed lines connect the mean values in 10 phase bins 
of the modulation. 



panel shows already hardly any loop structure wider than 
the scatter of the data. 



3.2 Changes in the global physical parameters 
during the Blazhko cycle 

In ISodor. Jurcsik fc Szeidll (|2009h we demonstrated that 
from good quality multicolour photometric data the physi- 
cal parameters of RRab stars can be determined with similar 
accuracy as with direct Baade-Wesselink analysis, using an 
inverse photometric method (IPM) . We applied this method 
to the BVIc light curves of MW Lyr at different phases of 
its Blazhko modulation succesfully, about 1 — 2% changes 
in the mean global physical para meters (L, T e f f , R) were 
detected during the Blazhko cycle (I Jurcsik et alj|2009l ). 

Before analysing the light curves of DM Cyg at different 
phases of its modulation, first we have to determine those 
parameters that definitely do not vary during the Blazhko 
cycle. These are the metallicity, the mass and the distance 
of the star and the dereddened standard magnitudes of the 
comparison stars i.e., the zero points of the magnitude scales 
used. 

The [Fe/H] of DM Cyg w as det ermined spectr oscop- 
ically by ISuntzeff et al.l (|l994f ) and iLavdenl (|l994h . The 
Fourier paramete rs of the mean V l i ght cu rve give [Fe/H] = 
—0.01 using the lJurcsik fc Kovacsi (|l996l . Eq. 3) formula. 
As the spectroscopic observations correspond to —0.16 
and 0.07 [Fe/H] values on the metallicity scale of the 
lJurcsik fc Kovacsi (|l996h formula, atmosphere models with 



[Fe/H ]= 0.0 and —0.1 compositions (|Castelli fc Kuruczl 
120031 ) are used in the light curve modelling of DM Cyg. 

The IPM finds mass values that are too large on evolu- 
tionary grounds in some cases if the mass is also allowed to 
vary by the fitting process. The fitting accuracy of the dif- 
ferent mass solutions differs, however, only marginally. The 
mean B,V,Ic light curves of DM Cyg are fitted with 5.1, 
4.7 and 4.0 mmag residual scatters when the mass value is 
also allowed to vary. In this case the best solution is found at 
m = 0.84SHs un . The r.m.s. of the fixed M = 0.55501^ mass 
solution B, V, Ic light curve fits are 5.2, 4.6 and 4.3 mmag, 
respectively. The average fitting accuracy in the three bands 
decreases only 0.1 mmag, from 4.7 to 4.6 mmag, when the 
mass is also allowed to be fitted. That is, the method is 
rather insensitive to the value of the mass. Therefore, in or- 
der to obtain reliable solutions, the IPM has been run with 
0.50, 0.55 and 0.60 9Jlsun fixed mass values for the mean 
light curves of DM Cyg. 

The global mean absolute physical parameters of DM 
Cyg averaged over both the pulsation and the Blazhko cycles 
and its distance are derived from the mean light curves using 
the IP method. The results are summarized in Table [3] as- 
suming different possible mass and [Fe/H] values. Note that 
the zero points of the colour and magnitude scales defined 
by the dereddened magnitudes of the comparison stars, are 
also allowed to vary i.e, they are also fitted at this step. To 
derive the distance, the apparent dereddened mean V mag- 
nitude of DM Cyg has to be known as the IPM gives the 
absolute mean V brightness of the star. To determine the 
apparent dereddened mean magnitude of DM Cyg, the stan- 
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Table 4. Mean physical parameters of DM Cyg derived from its mean light curves using the IP method. 
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fixed 
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Figure 12. Variations of the observed mean (left-hand panels) and derived parameters of DM Cyg during the Blazhko cycle. Circles and 
crosses denote the results when the IP method use Liu's template radial velocity curve, and when the V ra d curve is calculated from the 
Iq light curve, respectively. No significant, systematic difference between the results can be seen depending on the choice of the initial 
radial velocity templates. The results of the IP method for fixed 0.50, 0.55 and 0.60 OT/3Jt Slln mass values show the same variations in 
each parameters. Only the absolute mean values of the parameters are different, according to the data given in Table [4] These results 
prove that the detected changes of the physical parameters during the Blazhko cycle of DM Cyg are not sensitive to the settings of the 
IP method, and are independent of the uncertainties of the values of the mean global parameters. See further details in the text. 



dard V magnitude of the comparison stars and the interstel- 
lar absorption Ay should have to be known. Unfortunately, 
no standard magnitudes of the comparison stars have been 
published. We estimate the average of the dereddened V 
brightnesses of the two comparison stars to be 11.333 mag 
trans formed from their Tycho Bt, and Vt magnitudes <|ESAl 
1 19971 ) and using A v = 3.UE(B -V) = 0.558 interste llar ab- 
sorption value according to the ISchlegel et al.1 (|l998l ) maps. 
It has to be emphasized, however, that the uncertainties of 
the standard V magnitudes of the comparison stars and the 



interstellar absorption value affect only the distance esti- 
mate. Neither any other parameter nor their detected varia- 
tions depend on the choice of the dereddened V magnitudes 
of the comparison stars. 

The uncertainties of the estimates of the absolute phys- 
ical parameters listed in Table [4] correspond to the standard 
deviations of the results of running the IP code with 16 
different settings. These settings are the combinations of: 

i) B, V and Ic, or V, B — V and V — Ic data are used; 

ii) the V ra d curve is defined by Liu's template or it is cal- 



10 J. Jurcsik et al. 



culated from the Ic light curve; Hi) 2 different weights of 
the initial V ra d curve are used; iv) 2 different values of the 
AA(y ra d) I AVAmp ratio valid for Blazhko variables are ap- 
plied (see further d etails in ISodor. Jurcsik fc Szeidll |2009| ; 
I Jurcsik et aTll2009l ). However, as it can be seen from the 
data listed in Table [3] the true uncertainties of the mean 
global parameters are larger than the estimated errors of 
the solutions obtained for fixed mass and metallicity values. 
The true possible parameter space comprises the full range 
of the solutions for the possible mass and metallicity ranges. 
Nevertheless, we are focusing on the variation in the mean 
physical parameters and not on their absolute values, and it 
is shown in the rest of this section that the uncertainties of 
the mean global physical parameters have no effect on their 
variation during the Blazhko cycle. 

Though the amplitude of the modulation of DM Cyg is 
only about one tenth of that of MW Lyr, we have applied 
the IPM to the light curves in 10 different phases of the mod- 
ulation in order to decide whether there are any detectable 
changes in the mean global parameters of DM Cyg during 
its Blazhko period. The distance and the mass are fixed to 
their possible values given in Table [4] assuming [Fe/H]= 0.0 
metallicity. The zero points of the magnitude scales are also 
fixed correspondingly to the solution obtained for the mean 
light curves. 

Fig. [12] shows the results for the solutions using 9JI = 
0.509Jt Su n, d = 1179 pc; DJl = 0.55OT Stl n, d = 1215 pc and 
QTt = 0.609Jtsun, d = 1233 pc mass and distance values. In 
the left-hand panels, the variations in the observed mean 
quantities are plotted: total pulsation amplitude in V band, 
the variation of the pulsation period derived from the phase 
variation of the /o pulsation frequency and different pulsa- 
tion averages of the magnitudes and colours. The right-hand 
panels show the quantities derived from the IP method for 
three possible mass/distance combinations: the mean val- 
ues of the radius, the surface gravity, the absolute visual 
brightness averaged by magnitude and intensity units, the 
luminosity and the effective temperature. The different com- 
binations of the mass and distance give very similar results, 
the amplitudes and phases of the detected changes in the 
mean physical parameters during the Blazhko cycle hardly 
change, only their averages over the Blazhko cycle are dif- 
ferent, corresponding to their respective values given in Ta- 
ble[4] If atmosphere models with [Fe/H] =— 0.1 are used, the 
results on the variations of the physical parameters during 
the Blazhko cycle are not changed. 

In order to get a real estimate of the uncertainties of 
the derived quantities in different phases of the modulation, 
the IP code has been run with 16 different settings for each 
data set again. Note that the zero points of the colour and 
magnitude scales defined by the dereddened magnitudes of 
the comparison stars, are also fixed now. These magnitudes 
were, however, allowed to vary i.e., they were also fitted 
when the mean global parameters were determined. This 
explains why the scatter of the results for the 16 different 
setting of the IP method is much smaller in the different 
phases of the modulation than the scatter of the results for 
the mean light curves. 

The different intensity and magnitude averages of the 
V light curve and the colour curves show hardly any vari- 
ations, the amplitudes of their changes is only 0.002-0.004 
mag. The V, and B — V averages show some small system- 



atic variations of similar character like these averages in MW 
Lyr, but the V — I averages show only scatter. In spite of 
that the observed mean magnitudes and colours vary only 
slightly during the Blazhko cycle in DM Cyg, by the aid 
of the IP method 0.3% and 7 K systematic changes in the 
mean luminosity and temperature of the star could be de- 
tected during its Blazhko cycle. The phase relation of these 
variations are the opposite, similarly as it was detected in 
MW Lyr. Both stars are the most luminous and the coolest 
at around the largest amplitude phase of the modulation. 

The IP method calculates the mean physical parame- 
ters as their mathematical averages over the pulsation cycle. 
F($) =< F(<p) >, where $ and tp denote modulation and 
pulsation phases, respectively, and <> denotes averaging by 
ip for the whole pulsation cycle. For the luminosity, it means 
that the values plotted in Fig. 1 121 correspond to: 

L($) = < Hip) > = < 4naR(ip) 2 T eff {ip) A >, 
as the IP method satisfies the Stephan-Boltzmann law in 
eac h phase of the p ulsatio n (ip). 

ICarnev et all (|l992l ) showed from direct Baade- 
Wesselink analysis results that, even when the equilibrium 
luminosity and radius of pulsating variables equal with the 
arithmetic means of their variations over the pulsation cy- 
cle, the equilibrium temperature T eq is not the same as the 
mean temperature < T e // > of the star. The equilibrium 
temperature is defined then as 

T eq = (L eq /4naR 2 eq )^ 4 = (L($)/4^aJ?($) 2 ) 1 / 4 
= (< L{<p) > /47ia < R{ip) > 2 ) 1/4 . 
Consequently, the relation between the arithmetic means of 
the physical parameters shown in Fig. [12] deviates to some 
extent from the theoretical Stefan-Boltzmann law. 

The 0.13% variation in the mean radius is in good agree- 
mant with the 0.17% changes of the pulsation period (the 
pulsation equation requires AP/P ~ 3/2AR/R). The pul- 
sation period changes are determined from the phase differ- 
ences of the /o pulsation frequency in different phases of the 
modulation, but when deriving the radius variation the IP 
method does not utilize this information in any way. Con- 
sequently, the determined period and radius variations are 
completely independent quantities. The fact, that they show 
such a good agreement is a great support both to the IP 
method, and to our interpretation of the phase modulation 
as variations in the pulsation period. 

Comparing the res ults obtained for DM Cyg with that 
of MW Lyr (Fig. 14 in I Jurcsik et alj|2009h the similar be- 
haviour of the two stars are conspicuous. The only difference 
is in the detected amplitudes of the variations of the differ- 
ent parameters, each being about one tenth of their detected 
amplitudes in MW Lyr, in accordance with the different am- 
plitudes of the modulation of the two stars. 



4 THE PHOTOGRAPHIC AND 
PHOTOELECTRIC DATA 

The p hotographic light curve of DM Cyg was shown in lHurtal 
(2009). Although the Konkoly photographic data are very 
sparse, the observations covered 24 years with a 13-year gap, 
we tried to analyse this data set also searching for any sign 
of light curve modulation. No indication of light curve mod- 
ulation has been found in any of the two parts of the pho- 
tographic data, but there is some hint that modulation fre- 
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Figure 15. The O — C diagram of DM Cyg. Literature data are 
from the GEOS database (dots), the most deviant data are omit- 
ted. The Konkoly photographic, photoelectric and CCD normal 
maximum timings listed in Tabic [2] are plotted by squares, the 
time intervals of the Konkoly observations are indicated by hori- 
zontal lines in the top of the figure. The O — C indicates steady 
period increase during the ~ hundred years of the observations. 



Figure 13. Amplitude spectrum of the residual photographic 
light curve (after the removal of the pulsation signal) in the vicin- 
ity of the /o,2/o and 3/o pulsation frequencies. The residuals 
show signals at around /o + f m , 2/o + f m and 3/o + fm frequen- 
cies with 0.02 — 0.03 mag amplitude. The fo and f m frequencies 
detected in the photographic data are 2.3817564 and 0.0940, re- 
spectively. Bottom panel shows the spectral window of the pho- 
tographic observations. 
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Figure 14. Maximum brightness and phase of the B pg = 
11.5 mag brightness values on the rising branch (0.2 mag below 
maximum brightness) for the photographic observations phased 
with the supposed modulation period (10.636 days) are shown in 
the left-hand and right-hand panels, respectively. Data from the 
two parts of the photographic observations are shown with dif- 
ferent symbols. According to these plots, there is some indication 
that the maximum brightness and the phase of the rising branch 
of DM Cyg varied with small amplitude during the time of the 
photographic observations as well. 



quency components appear at fc/o + fm, (k = 1, 2, 3) posi- 
tions in the residual spectrum of the complete photographic 
data set (see Fig. I13[l . The highest modulation peaks are not 
exactly at the same separations for fo + f m , 2/o + f m and 
3/o + f m - The best solution which gives the smallest resid- 
ual for the photographic data can be gained with 2.3817564 
c/d (0.41985822 d) and 0.0940 c/d (10.636 d) pulsation and 
modulation frequencies, respectively. 

Fig. [14] plots the maximum brightness and the phase 



phased with the supposed 10.636 d modulation period. The 
phase of a given magnitude on the rising branch can be de- 
termined with higher accuracy than the phase of the max- 
imum, therefore we use this quantity to measure the am- 
plitude of the phase modulation component in the photo- 
graphic data. Data belonging to the two parts of the obser- 
vations are denoted by different symbols in Fig. Q3] Accord- 
ing to these plots the amplitudes of the amplitude and phase 
modulations were not larger at the time of the photographic 
observations than today. 

The photographic data were also analysed taking into 
account the period decrease of the pulsation that took place 
during the time interval of the observations, but the results 
were only marginally different from that obtained from the 
original photographic data. 

The photoelectric observations comprise data only from 
four nights which does not allow to check the light curve 
modulation. The photoelectric data were therefore used only 
to determine one normal maximum timing listed in Table [2] 

The normal maximum timings of the Konkoly obser- 
vations given in Table [2] together with literature data col- 
lected in the GEO S database are plotted in Fig. 1151 Based 
on the GEOS data lLe Borgne et all (|2007T ) determined p = 
0.091d Myr -1 period increase rate for DM Cyg. The addi- 
tion of the Konkoly data does not modify this result. 

From the analyses of the CCD and photographic data of 
DM Cyg we conclude that the period changes of the pulsa- 
tion and modulation have the opposite direction, when the 
pulsation period of DM Cyg was shorter by 0.000005 d then 
the modulation period was 0.066 d longer than today. These 
values correspond to ^ = -13200 or ^pf/^ = -523 pe- 
riod change rates. However, these data have to be taken with 
caution partially because the uncertainties of the results de- 
rived from the photographic data and also because pulsation 
and modulation period changes of those Blazhko variables 
where the periods could be determined for several epochs 
show that there is no strict relation between the two periods, 
for some time intervals even the sign of the period change 
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ratio can change (e.g., in RV UMa lHurta et al" . 2008]). There- 
fore, the period change ratios defined by only two-epoch data 
should be missleading in some cases. We also have to note, 
however, that the complex period change behaviour of the 
modulation is usually connected to complex changes of the 
pulsation period and of the pr operties of the modulation, too 
(e.g., in RR Gem and XZ Cve lSodor, Szeidl fe Jurcsikll2007l ; 
lLaCluvze et alj|2004h . Probably, the steady period change 
of the pulsation of DM Cyg is accompanied with steady pe- 
riod change of its modulation and with the stability of its 
modulation properties. 



5 SUMMARY 

In the present paper the light curve modulation of DM Cyg, 
a fundamental mode RR Lyrae variable has been investi- 
gated. Most of the results are very similar to those found 
in the analysis of MW Lyrae, a large modulation amplitude 
Blazhko variable, but on a significantly smaller scale. The 
0.07 mag amplitude of maximum brightness variation of DM 
Cyg is only 10% of the modulation amplitude of MW Lyr, 
accordingly, all the changes detected in DM Cyg during its 
Blazhko period are about 10% of the detected changes in 
MW Lyr. 

The phase and amplitude relations of the amplitude 
and phase modulation components are similar in DM Cyg 
and MW Lyr. The variations of the mean magnitudes and 
colours in DM Cyg where they can be detected are also sim- 
ilar to their counterparts in the MW Lyr data. Therefore, it 
is not surprising that the derived changes in the mean global 
physical parameters of DM Cyg resemble the changes found 
in the mean parameters of MW Lyr only on a much reduced 
scale. 

The only notable difference between the character of 
the modulations of the two stars is found when the changes 
of the Fourier parameters of the light curves during the 
Blazhko cycle is investigated. While in the case of MW Lyr 
the phases of 2/o, and 3/o change in line with the phase vari- 
ation of /o (i.e, there is no significant changes detected in the 
f(fki) phase differences), when analysing the light curves of 
DM Cyg 90° phase differences between the phase variation 
of the /o pulsation frequency and the phase variations of 
the higher order pulsation components are found. In spite of 
this, the data transformation that corrects the times of the 
observations taking into account the phase variation of the 
/o pulsation frequency which completely separated the am- 
plitude and phase modulation components of the light curve 
modulation of MW Lyr seems to work also on the light curve 
of DM Cyg. 

The modulation behaviour of other Blazhko variables 
can differ, however, more significantly from the modulations 
of DM Cyg and MW Lyr. There are Blazhko stars, where the 
phase relation of the amplitude and phase modulation com- 
ponents contrast with the phase relations of DM Cyg and 
MW Lyr. Also, the detected changes in the mean colour s 
may vary differently e.g, in SS Cnc (| Jurcsik et al.l l2006h . 
Therefore, we cannot draw major conclusions from the sim- 
ilarity of the results obtained for DM Cyg and MW Lyr be- 
fore analysing other Blazhko variables of different character 
on a similar way. 
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